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Abstract 

We propose that the CP violating phase in the CKM mixing matrix is identical to the CP phases 
responsible for the spontaneous CP violation in the Higgs potential. A specific multi-Higgs model 
with Peccei-Quinn (PQ) symmetry is constructed to realize this idea. The CP violating phase 
does not vanish when all Higgs masses become large. There are flavor changing neutral current 
(FCNC) interactions mediated by neutral Higgs bosons at the tree level. However, unlike general 
multi-Higgs models, the FCNC Yukawa couplings are fixed in terms of the quark masses and CKM 
mixing angles. Implications for meson-anti-meson mixing, including recent data on D — D mixing, 
and neutron electric dipole moment (EDM) are studied. We find that the neutral Higgs boson 
masses can be at the order of one hundred GeV. The neutron EDM can be close to the present 
experimental upper bound. 



1 



I. INTRODUCTION 



The origin of CP violation is one of the outstanding problems of modern particle physics. 
There have been several experimental measurements of CP violation ll. All of them are 

d r, 

consistent with the Cabbibo-Kobayashi-Maskawa (CKM) model [2j, |3j, where the source 
of CP violation comes from the phase [3|] 5km in the CKM mixing matrix for quarks. A 
successful model of CP violation at the leading order should have the successful features of 
the CKM model. It is important to understand the origin of CP violation. An interesting 
proposal due to T.-D. Lee was that CP is spontaneously violated Id] . The popular Weinberg 
model [5I of spontaneous CP violation model has problems Q, 7] with data and has been 
decisively ruled out by CP violating measurement in B decays Spontaneous CP violation 
in the Left-Right models has also been ruled out for the same reason^]. In this work we 
restore the idea that CP is broken spontaneously and the phase 5km is the same as the 
phase 5 spon that causes spontaneous CP violation in the Higgs potential. We construct 
specific models to realize this idea. The main difference of our models lies in how the CP 



violating phase in the CKM matrix is identified 9J. 

Let us start by describing how a connection between 5km and 5 spon can be made. It is 
well known that to have spontaneous CP violation, one needs two or more Higgs doublets 
0j. Consider the following Yukawa couplings with multi-Higgs doublets, 



L Y = Qi(r«i0i + ^u2<t>2)U R + QiXd^dDa + h.c. 



(1) 



where Ql, Ur and D R are the left-handed doublet, right-handed up and right-handed down 
quarks, respectively. Generation indices are suppressed, (pd = —i&2<f>d an d </><2 m ay be one of 
the 0i 2 or another doublet Higgs field. The Yukawa couplings r n i 5 «2,d must be real if CP is 
only violated spontaneously. 

The Higgs doublets when expressed in terms of the component fields and their vacuum 
expectation values (vev) V{ are given by 

^ — = 1 r . -i- 1 i -i- / . -t . i 1 

(2) 



-^(vi + Ri + iAi 



V 



The quark mass terms in the Lagrangian are 



-Ur 



M ul e^ + M u2 e^ U R - D L M d e~^D R + h.c. 



J82 



(3) 
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where M ui = -T ui Vi/\/2. 

The phases 9\ and 9d can be absorbed by redefining the fields Ur and Dr. However, the 
phase difference 5 = 62 — 61 cannot be removed and it depends on the Higgs potential. A 
non-zero 5 indicates spontaneous CP violation, 5 = S spon . Without loss of generality, we 
work in the basis where Dl, Dr are already in their mass eigenstates. In this basis the 
down quark mass matrix is diagonalized, which will be indicated by M^. In general the 
up quark mass matrix M u = M u \ + e tS M u2 is not diagonal. Diagonalizing M u produces the 
CKM mixing matrix. One can write M u = VckmM u V r . Here Vckm is the CKM matrix 
and Vr is an unknown unitary matrix. A direct identification of the phase 5 spon with the 
phase 8km m the CKM matrix is not possible in general at this level. There are, however, 
classes of mass matrices which allow such a connection. A simple example is provided by 
setting Vr to be the unit matrix. With this condition, M u = Vq KM M u . One then needs to 
show that Vq KM can be written as 



Vi KM = (M ul + e iS M u2 )M-\ 



(4) 



Expressing the CKM matrix in this form is very suggestive. If Vckm (or Vckm) can always 
be written as a sum of two terms with a relative phase, then the phase in the CKM matrix 
can be identified with the phase 5. 

We now demonstrate that it is the case by using the Particle Data Group (PDG) 
parametrization as an example. To get as close as to the form in Eq. (jljh we write the 
PDG CKM matrix as fl 

" " ' 12' ].:!' -M2 1 I.:-!' -M:-S ^ 

Vckm 





1 
1 



-S12C23 - Ci2S23Sl 3 e J<513 C12C23 - Sl2S23Sl 3 e 2 ' 513 



Sl3 

S23C13 



\ S12S23 - Cl 2 C23Si 3 e 



j<5l3 



(5) 



-C12S23 - Si 2 C23Si 3 e'" 1J C23C13 j 



where s,- 



sin 9ij and Cij = cos 6{j . 



Absorbing the left matrix into the definition of Ul field, we have 



\ 



( 



M„ , M, 



u2 



C12C13 — Ci2S23 s 13 — C12C23S13 
S12C13 — S12S23S13 — S12C23S13 





\ 



(6) 



I — S12C23 S 12 S 2 3 

C12C23 — C12S23 

V S13 S23C13 C23C13 J 

and 5 = —S^. We therefore find that it is possible to identify the CKM phase with that 
resulting from spontaneous CP violation. Note that as long as the phase 5 is not zero, CP 



violation will show up in the charged currents mediated by W exchange. The effects do not 
disappear even when Higgs boson masses are all set to be much higher than the W scale. 
Furthermore, M 12 are fixed in terms of the CKM matrix elements and the quark masses, as 
opposed to being arbitrary in general multi-Higgs models. 

We comment that the solution is not unique even when Vr is set to be the unit matrix. 
To see this, one can take another parametrization for the CKM matrix, such as the original 
Kobayashi-Maskawa (KM) matrix [3|. More physical requirements are needed to uniquely 
determine the connection. The phenomenological consequences will therefore be different. 
We will come back to this when we look at phenomenology of models. The key point we 
want to establish is that there are solutions where the phase in the CKM matrix can be 
identified with the phase causing spontaneous CP violation in the Higgs potential. 

The mass matrices M u \ and M u2 can be written in a parametrization independent way 
in terms of the eigen-mass matrix M u , the CKM matrix, and the phase 5, 

Mm = V ] CKM M U - ^-I m {vi KM )M u , 

M u2 = -±-Im(V^ KM )M u . (7) 
sin o 

Alternatively, a model can be constructed with two Higgs doublets couple to the down 
sector and one Higgs doublet couples to the up sector to have, 

L Y = Q L Tu<t>uU R + Q L (T dl 4>i + T d2 4> 2 )D R + h.c. . (8) 
In this case M di = —T di Vi/\/2, and 

M dl = V CKM M d + ^—Im(V C KM)M d , 
smo 

M d2 = - — Im{V C KM)M d . (9) 
sin o 

We denote the above two possibilities as Model a) with two Higgs doublets coupled to 
the up sector, and Model b) with two Higgs doublets coupled to the down sector. 

II. MODEL BUILDING 



A common problem 
term will be generated 



or models with spontaneous CP violation is that a strong QCD 9 



. Constraint from neutron dipole moment measurement will rule 
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out spontaneous CP violation as the sole source if there is no mechanism to make sure that 
the 9 term is small enough if not zero. The models mentioned above face the same problem. 



We therefore supplement the model with a Peccei-Quinn (PQ) symmetry 10] to ensure a 
small 9. 

To have spontaneous CP violation and also PQ symmetry simultaneously, more than two 
Higgs doublets are needed Jill]. For our purpose we find that in order to have spontaneous 
CP violation with PQ symmetry at least three Higgs doublets <pi = e tdi Hi and one complex 
Higgs singlet S = e t9s S = e l6s (v s + R s + iA s ) /y/2 are required. The Higgs singlet with a 
large vacuum expectation value (vev) renders the axion from PQ symmetry breaking to be 
invisible [l^ . thus satisfying experimental constraints on axion couplings to fermions. 



12|. 



We will henceforth work with models with an invisible axion 
The PQ charges for Models a) and b) are as follows 

Model a) Q L : , U R : -1 , D R : -1 , 0i, 2 : +1 , d = 3 : -1; 

Model b) Q L :0, U R : +1 , D R :+1, 0i )2 : +1 , <j> u = fa:-l. (10) 

In both cases, S has PQ charge +2. For leptons, the PQ charges can have different assign- 
ments. For example: L L : , e R : — 1 or L L : , e R : +1. 

For both models a) and b), the Higgs potentials have the same form which is given by 

V = -m\H\H 1 - m\H\H 2 - m\E\E z - m 2 l2 (HlH 2 e i{e2 - dl) + h.c.) - m 2 s S ] S 
+ X^HlH,) 2 + \ 2 {HtH 2 f + X t (HlH 3 ) 2 + \ S (S^S) 2 

+ A 3 (# 1 ^i)(# 2 t i/ 2 ) + A , 3 (iy 1 t F 1 )(iy 3 t F 3 ) + \'^hIh 2 )(hIh 3 ) 

+ \i(HlH 2 )(H 2 Hi) + ^{HlHs){HlH x ) + Xl(H j 2 H 3 )(HlH 2 ) 

+ h 5 ((HlH 2 ) 2 e 12 ^-^ + h.c.) + \ 6 (HlH 1 )(HlH 2 e>fr- < » + h.c.) 

+ \ 7 (HlH 2 )(H\H 2 e^~^ + h.c.) + \ % (HlH z )(HlH 2 c^- d ^ + h.c.) 

+ hH\H x S^S + f 2 H ] 2 H 2 S ] S + f 3 HlH 3 £fiS + d 1% {H\H 2 e i{ - ^ + HlH^-^-^S^S 

+ f 13 {HtH 3 Se^ 3+e ^ + h.c.) + f 23 (HlH 3 Sc t{e * +e °- e2) + h.c.) . (11) 

Only two phases occur in the above expression, which we choose to be 8 = 9 2 — 9\ and 
5 S = 9 3 + 9 S — 9 2 . The phase 9 3 + 9 s — 9\ can be written as 5 + 5 S . Differentiating with respect 
to S s to get one of the conditions for minimization of the potential, we get 

ii3^3^ sin(5 s + 5) + f 23 v 2 v 3 v s sin S s = . (12) 



We see that 5 and S. are related with 



tan 5, 



/i 3 ui sin 5 



(13) 



f23V 2 + /i 3 f 1 COS 5 ' 

Therefore, S is the only independent phase in the Higgs potential. A non-zero sin S is the 
source of spontaneous CP violation and also the only source of CP violation in the model. 

In this model the Goldstone fields h w and h z that are "eaten" by W and Z, and the axion 
field are given by 



h w = -(^i^i + v 2 h 2 + v 3 h 3 ) , 

h z = -{viAx + v 2 A 2 + v 3 A 3 ) , 
v 

a = (-vtvjA! - v 2 v\A 2 + v\ 2 v 3 A 3 - v 2 v s A s )/N a , 



(14) 



where v 2 = v \ + v 2 + v\ and N 2 = (vj 2 v 3 v 2 + v A v 2 s ) with v 
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We remove h w and h z in the Yukawa interaction by making the following changes of basis, 



ua 




v 2 /v 12 


-viv 3 v s /N A 


vi/v 


-viv 2 /N a ^ 


(aA 


A 2 




-v 1 /v 12 


-v 2 v 3 v s /N A 


v 2 /v 


-v 2 v 2 /N a 


a 2 


A 3 







vf 2 v s /N A 


v 3 /v 


v 2 2 v 3 /N a 


h z 


[a s j 




V 


vf 2 v 3 /N A 





-v 2 v s /N a J 


V a ) 



h 2 
\h 3 / 



( 



V 2 /V 12 V&3/VV12 vi/v 



\ 



f Hr \ 



(15) 



-V1/V12 v 2 v 3 /vv 12 v 2 /v H 2 
\ -v 12 /v v 3 /v j \ h w ) 

where N\ = v 2 2 (v 2 2 v 2 + v 2 s v 2 ). a 1:2 and H{ 2 are the physical degrees of freedom for the 
Higgs fields. With the same rotation as that for the neutral pseudoscalar, the neutral scalar 
Higgs fields (R u R 2 , R 3 , R S ) T become (if°, H 2 , H®, H^) T . Since the invisible axion scale i>, 
is much larger than the electroweak scale, to a very good approximation, N a = v 2 v s and 
N A = v 12 vv s . 

In the rotated basis described above, we have the Yukawa interactions for physical Higgs 
degrees of freedom as the following 

AS 



4 a) = U L [M U - 



^--(M u - VcKM Im{Vl KM )M u ^-)^- 
v 12 v 2 smd v\v 2 



\U R {H» + ia\) 



v 12 v 



- D L M d D R [^(H° - ia 2 ) + -H° 3 

V 3 V V 



'12 



{Hi - ia)\ 



+ V2D L [V£ KM M„ 



iS 



Vl2 

('2('j 2 smS / Viv 2 

;Vl2 



~ (V6kmM u - Im{V J CKM )M u - 



V2^-D L V£ KM M U U R H; - V2^U L V C KMM d D R Ht + h.c. , 
Vi 2 V ~ vv 3 



D L [M a 



V12V2 



~(M d 



V CKM 



Im{V C KM)M d ^-)^\D R {Hl - m?) 



+ D L M d D R [^{H Q 2 - ia 2 ) - -H° + 

V 12 V V 

1 



V 2 V , 



sin 5 v \v 2 
[Hi - ia)\ 



U L M u U R p^{H° 2 + ia 2 ) + -H° 3 

V 3 V V 



U 12 

V 2 Ve 



'M + ia)] 



-iS 



- V2U L [V CKM M d ^- - {VcKMMd + Im{V C KM)M d ^—)^-]D R Hl 



V2V12 



sin 5 V1V2 



+ V2^U L V C KMM d D R H+ + V2^D l V^ km M u UrH 2 + h.c. 



V12 



(16) 



vi 2 v ~ vv 3 

Note that the couplings of a and H% to quarks are suppressed by l/v s , and that only the 
exchange of and a\ can induce tree level FCNC interactions. The FCNC couplings are 
proportional to VcKMlm{VcKM)M u and VQ KM Im(VcKM)M d for models a) and b), respec- 
tively. 

We have mentioned before that the identification of the phase 5 with that in the CKM 
matrix does not uniquely determine the full Yukawa coupling. Here we give two often used 
parameterizations, the PDG CKM matrix and the original KM matrix with the CP violating 
phase indicated by 5 km, to illustrate the details. In the two cases under consideration, the 
phase 5 are identified with — 5± 3 and — 5 km, respectively. The differences will show up in the 
FCNC of neutral Higgs coupling to quarks which are proportional to the following quantities, 



PDG : VcKMlm(y£ KM )M v = - sin 5 13 e^ :i 



'13 



"S23 s 13 c 13 



— S23S13C13 — C23S13C13 



„2 2 
* 23*13 



\ — c 23 s 13 c 13 S 23 C 23 S 13 

' c\ 2 S12C12 ^ 



V ] CKM Im{V CKM )M d = sm5 13 e- iS ™ 



Sl2 c 12 



/ 



'12 








KM : V C KMlm{y^ KM )M u = - sin 5 K m& 



i5 K M 








-S2C2 



\0 -s 2 c 2 c\ j 



523 C 23 S 13 

c 2 2 
c 23*13 



M„ 



7 



VcKM Im { v CKM)M d = sm5 KM e 



-i&KM 





s| -S3C3 
^0 -S3C3 c\ 



M, 



d ■ 



(17) 



III. MESON AND ANTI-MESON MIXING AND NEUTRON EDM 

In this section we study some implications for meson and anti-meson mixing and neutron 
electric dipole moment. 



A. Meson and Anti-meson Mixing 



Meson and anti-meson mixing has been observed previously in K° 



11] 



and in D° — D° very recently 14|. In the models considered in the previous section, besides 
the standard "box" diagram contributions to the mixing due to W exchange, there are also 
tree level contributions due to the FCNC interactions of H® and a\. 

The interaction Lagrangian for Hi and a k with quarks have the following form for both 
models a) and b), 

L = qMj + blj^QjH? + iqi{c% + d\f ih )q i a k . (18) 

For the meson and anti-meson mixing, only the FCNC interaction terms of H\ and a\ 
contribute. We can write a 1 = d 1 = a and b l = c 1 = j3, with a = (A + A^)/2 and 
fi= (A - v4 f )/2, and A given by: 



For a), A 



t e V 12 

VcKMlm ( Vckm) M u 

sm V1V2 



For b), A = -V c 



' Im(V C KM)M d e n2 



(19) 



CKM " "V ' l/^ju; u . j 

suit) t>if 2 

Using the definition < 01^7^75^- >= ifppp/ y/2mp and the equation of motion q^qj = 
(j>i ~ PjYqtl 'ixlhQ.3 1 '( m i + m j) with p p = pj — Pi, we obtain the matrix element for P — P 
mixing in the vacuum saturation approximation as 

f 2 P m] 



M- 



12 



in 



Hi 



mi 



(7,2 _ _f a 2. , 6 2 ^ JP'"P 

12 1 ^ + ^ jj K + m,) 2 

3 ?. _ J_C a 2. , fe 2 ^ /p m P 



+ j^(b 2 ij- a2 ij)fp m P 



+ ^(Oy -kj)fpmp 



(20) 



m 2 Hl m ai 6 (mj + mj) 2 



where m 2 Hiai parameterizes the mixing between a% and Hi, that is determined from the Higgs 
potential V = m 2 Hiai Hiai + .... Since rn 2 H iai involves new parameters, it can be made small 
enough to avoid any conflict with data. We will neglect their contributions when discussing 
meson and anti- meson mixing. We will come back to this when discussing neutron EDM. 

It is obvious from the structure of the flavor changing coupling in Eq. (117j) that the flavor 
changing structure for the PDG and KM cases are different. For PDG case, in model a) 
there is tree level contribution from neutral Higgs exchange to D° — D° mixing, but no 
contribution to K°, B® and B° s mixing. In model b), there is only non-zero contribution to 
K° — K° mixing at the tree level. For the KM case, there is no tree level contribution to 
meson mixing in model a). For model b), there is only non-zero contribution to B® mixing. 

In our numerical analysis, we will use the following values for the relevant parameters. 
For the CKM matrix elements, we take the PDG central values with [1]: s 12 = 0.227, 
s 23 = 0.042, sia = 0.004 and sin<5 13 = 0.84 (equivalently Si = 0.227, s 2 = 0.0358, 
S3 = 0.0176 and sin<5 = 0.97 for the KM parameterization). For the quark masses, we 
take [16J m M (l GeV) = 5 MeV,m d (l GeV) = 10 MeV,m s (l GeV) = 187 MeV,m c (m c ) 
1.30 GeV, m^mt) = 4.34 GeV, m t = 174 GeV. For the meson decay constants, we take 
f K = 156 MeV, f D = 201 MeV, and f Bs = 260 MeV. 



Models with PDG parameterization 

We consider the models with PDG parameterization first. 

Model a): In this case there is mixing only in D° — D° system. Combining the BaBar and 



Be 



le [y| results the 68% C.L. range for x = Am/To is determined to be (5.5 ± 2.2) x 10 
15( | . Theoretically the elements in A for this particular case are Ayi = — S23S13C13 ^7^7 and 
A 2 i = -S23Si3Ci3 !Z ^f , which implies that a X2 ~ 612 ~ -s 2 3Si3Ci 3 f^f. We obtain 

5 2 2 2 ,vi 2 m c 2 f%m D m D 2 1 1 



V\V 2 



12 ~23^i3v ViV2 j Td y mc + m j ^ m 2j 

7 - 5 x 10" % ■ L 2 2 (—2 V)(100 GeV) 4 . (21) 



(sin2/?) 2 wf 2 m 2 Hl 



m 



where tan/3 is defined to be V\jv 2 . 

It is well known that the SM short distance contribution to the D — D mixing is 
small. Long distance contributions can be much larger, but they suffer from considerable 
uncertainty. New physics may contribute significantly |15 ]. It is tempting to see if the 
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new contribution in this model can account for the full measured value. If the effective 
neutral Higgs mass = l/(l/m 2 Hl — l/m? ai ) is of order 100 GeV, one would require 

sin 2 2/3t> 2 2 ~ (12) 2 GeV 2 . Since are related to the top quark mass, with the assumption 
that the top quark Yukawa coupling y t < 1, one of them should be large, ~ 240 GeV. 
Saturating the experimental central value for x, we would have sin(2/3) ~ 0.05 implying 
v\/v2 or V2/V1 to be of the order of 1/40. If all vevs are the same order of magnitude, the 
new contribution does not produce large enough x to saturate the measured value. 

Model b): In this case there is mixing only in K° — K° system. We have 



^ = 4.4xl0- 12 1 (4-1 
m K sin 2pv{ 2 m 



4 - 4 >< 1Q - 12 ^2o^. 2 (— - -r)( 100 GeV ) • ( 22 ) 

Hi ai 



This is to be compared with the experimental data Amx/mx = 7.0 x 10~ 15 . It puts strong 
constraints on the scalar masses, i.e., the Higgs particles must be at least TeV scale to 
suppress the value if ai and Hi are not degenerate in mass. 

Models with KM parameterization 



We now come to models with the original KM parameterization. In this case, there is no 
meson and anti- meson mixing in Model a). 

Model b): There is mixing only in B s — B s system. We have 

^ x io - i2 s?W4 - ^ )(io ° GeV)1 (23) 

Experimental value Am^ = 17.4 ps -1 implies Amgjm^ = 2.1 x 10~ 12 . It has been 



shown in Ref. 18j that the New Physics contribution to Am^ s can be up to 10%. To obtain 
the lowest Higgs boson mass, we maximize sin 2/3=1 which requires v\ = v 2. Taking f 1,2,3 to 
be all equal, the Higgs boson mass can be as low as 300 GeV. With smaller v 1,2 or non-equal 
v 1 and f 2 would lead to more stringent bound on Higgs mass. 

B. The neutron EDM 

The neutron EDM can also provide much information on the model parameters. The 



standard model predicts a very small 19j d n (< 10 31 e cm). The present experimental 
upper bound on neutron EDM d n is very tight [1]: \d n \ < 0.63 x 10~ 25 e cm. We now 
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study whether neutron EDM can reach its present bound after imposing the constraints 
from meson and anti-meson mixing discussed in the previous section. 

In the models we are studying, the quark EDMs will be generated at loop levels due to 
mixing between Oj and H^. From Higgs potential given earlier, we find the mixing parameters, 

m H iai = [(^6 - h)v\v 2 - \{v\ - v\) cos 5] sin 5 , 



m ka 2 



f 13 sin(5 + 8 s )vv s 



V2v 2 

m Lai — n [-2\ 5 v 1 v 3 vlsm25 + 2(-\ 6 vf - X 7 vj + (X 8 + d 12 )vf 2 )v 2 v 3 smS 

+V2f 13 v 2 v s sin(5 + 5 S )} , 
m 2 H3ai = ^-[2\ 5 viv 2 cos(6'i - 9 2 ) + X 6 v 2 + X 7 v% + X 8 vf\ sin 5. (24) 

Note that all the parameters above are zero if sin 5 = 0. 

The one loop contributions to the neutron EDM are suppressed for the usual reason 
of being proportional to light quarks masses to the third power for diagram in which the 
internal quark is the same as the external quark. In model a) with PDG parameterization, 
there is a potentially large contribution when there is a top quark in the loop. However, 
the couplings to top are proportional to Si 3 , therefore the contribution to neutron EDM is 
much smaller than the present upper bound. We will not discuss them further. 

It is well known that exchange of Higgs at the two loop level may be more important 
than the one loop contribution, through the quark EDM OjJ {^(J , quark color EDM O q 
and the gluon color EDM |2lj] defined as 



20|, 



V q = -^iqa^F^q , O c q = -^qa^G^q , O c g = —C f ahc G%G b m G c va , 



where is the photon field strength, G^ v is the gluon field strength and G^ u = -^e^^G 01 ^ . 

In the valence quark model, the quark EDM and color EDM contributions to the neutron 
EDM d n are given by Q 



d l = Vd ^d d - ^d u 



, d c n = erj f 

A 



4 2 



A 



(26) 



where U V, = ^ ^ ^ * 0.166 and „ = 

(m) Um (Sfei) 11 ' 25 (w) 11/27 m " °' 0117 are the <**> running faetors from 
scale mz to the hadron scale A. 
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A naive dimensional analysis (NDA) estimate gives the gluon color EDM contribution to 
the neutron EDM as the following 

d n « (27) 



where M = 4-7T = 1190 MeV is the scale of chiral symmetry breaking. The QCD running 
factor is H e = (-i-i)- 54/23 (-i-l)- 54/25 (^AL)- 54/27 « 1.2 x 10" 4 . 

! 'J s V 4?r / \ct B (m t ) J \a B (m b ) ) \a s (m c ) J 

The two loop contribution to d g , / 9 and C are given by 

dq = ^l mqG{q) , /g = ^G(g) , C = lH(g) , (28) 
where is the charge of q quark and 
G(q) 



mii mi mi, mi H 

til a k til O-k 



H{g) = (h(^f) - h(^-))ImZ% , (29) 
mi, mi 

til a k 

where ImZij is defined through ImZj* = 2a l ii djjXik/(miTrij) with a l ,d h defined by Eq. (fl8"j) 
and Aj& = m 'ji iak /(' m % l ~ m a fe ) i s a mixing factor depending on the neutral Higgs bosons 
exchanged in the loop. 

The functions f(z), g(z) and h(z) are given by 

z f 1 , l-2x(l-x) n x(l-x) 
= - / ax — -; - In 



2 Jo x(l — x) — z z 

9{z) = o / ^— r ln^ ^ , 30 

2 Jo x(l — xj — z z 

h(z) = - f 1 dx tdu A3(1 ~ X) 2 . (31) 

V; 2 Jo Jo [zx(l -ux) + (1 -u){l-x)f V ; 

Numerically we find that functions (/, g, h) change slowly from (0.5, 1, 0.1) to (0.2, 0.2, 0.03) 
when Higgs masses are increased from 100 GeV to 1 TeV. 



Models with PDG parameterization 

Model a): The 2- loop contributions to the neutron EDM due to the Higgs bosons exchange 
in the loop are proportional to the mixing factor Xik{f, g, h) . We take these factors to be 
approximately equal to estimate the contributions from different Higgs exchange. 

If using the parameters which produce D mixing, i.e., tan/5 = 40, V\2 ~ 240 GeV and 
i>3 ~ 10 GeV and Higgs around 100 GeV are used, we find that the dominant contribution 
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is from H 3 , a\ exchange, 

2 

d n « -1.5 X 1Q-25 m g 3 ax g cm _ ,3^ 

If all vevs are of the same order, i.e., taking v\ — v 2 — V3 with Higgs mass of order 100 
GeV, we have 

4~8x 10~ 26 9 H * ai - e cm . (33) 

m H 3 ~ m l x 

Model b): In this case H 1 ,ai do not couple to it, so the two loop contribution to quark 
EDM and quark and gluon color EDM from the Hi, a\ are small. 

The contributions to neutron EDM are about the same from the H x ,a 2 and H 23 ,ai ex- 
change, with different mixing factors. Explicitly as an example, for the case Hi, a 2 exchange 
with the Higgs mass taken to be 1 TeV, as high as allowed by K° — K° mixing, we have 

2 

1Q -26 m Hia 2 e cm _ ^ 
m k ~ m a 2 

If rn 2 Hia2 is not too much smaller than m 2 Hi a2 , the neutron EDM can be close to the upper 
bound. 

Models with KM parameterization 

Model a): In this case there are no constraints from meson mixing, the Higgs mass can 
be low. If all vevs are of the same order, i.e. taking v\ = i> 2 = 1*3 with Higgs mass of order 
100 GeV, we have the main contribution come from Hi, a 2 exchange, 

2 

d n W 5 X 1Q-26 m Hia 2 Q cm _ ^ 

m 2 Hi - m 2 a2 

Model b): Similar to the case for Model b) as in the PDG parameterization case, the 
contributions from the H\,a\ exchange are small. Taking the vevs to be same order and 
Higgs mass to be of the order of 100 GeV, we find the contributions from H\,a 2 exchange 
and H 2t 3,ai exchange are comparable. For the case H\,a 2 exchange, the contribution is 
given by 

2 

4~5x 10~ 26 9 mglQ2 - e cm . (36) 

If one takes the Higgs mass to be 300 GeV as that from B s — B s mixing, the neutron EDM 
will be smaller. 
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IV. DISCUSSIONS AND CONCLUSIONS 



In our previous discussions, we have not considered Yukawa coupling for the lepton sector. 
An analogous study can be carried out. If one introduces right handed neutrinos, see-saw 
mechanism can be applied to generate small neutrino masses if the axion scale v s is identified 
with the see-saw scale. We briefly discuss two classes of models parallel to the quark sector 
before conclusion. 

Model a): The PQ charges for lepton doublet Ll, electron e R and neutrino v R are assigned 
as: L L (0), e R (—l) and v R {— 1). The Yukawa couplings are then 

L = L L {Y X H X + Y 2 H 2 e i5 )u R + L L Y 3 H 3 e R + v c R Y s Se l{s+ ^v R + h.c. (37) 

In this case the mass matrices in L m = —e L M e e R — v L M D v R — (l/2)u < ^M R i , R can be written 

as 

M l = -^=Y 3 v 3} M D = -^=(Y lVl + Y 2 v 2 e l5 ), M R = —\f2Y s v s e ilyS+Ss \ (38) 

The charged current mixing matrix in the lepton sector, the Pontecove-Maki-Nakagawa- 
Sakata (PMNS) matrix 24], Vpmns similar to the Vckm matrix is given by Vpmns — ViXl 



i 

where V£ and V£ are defined by: M e = v£M e V% and M u = -M D M R l Ml = Vj?M v V£* 
with M e and M v the charge lepton and light neutrino eigen-mass matrices. 

Model b): The PQ charges for lepton doublet L L , electron e R and neutrino v R are assigned 
as: Ll(0), e R (+l) and u R (+l). The Yukawa couplings are 

L = L L Y 3 H 3 u R + L L (Y\H\ + Y 2 H 2 e^ 5 )e R + ^Y s S^-^h R + h.c. , (39) 

and 

M l = -^=(Y 1 v 1 + Y 2 v 2 e^ 5 ), M D = --^=Y 3 v 3 , M R = -V2Y s v s e^ s+ ^ . (40) 

From the above we see that, in general there are CP violation in the mixing matrix 
Vpmns-i an d the source is the same as that in the Higgs potential. But the identification of 
the phase 5 with the phase in the Vpmns becomes more complicated due to the appearance 
of Mr. The related details will be discussed elsewhere. 

We have proposed that the CP violating phase in the CKM mixing matrix to be the 
same as that causing spontaneous CP violation in the Higgs potential. Specific multi- 
Higgs doublet models have been constructed to realize this idea. There are flavor changing 
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neutral current mediated by neutral Higgs bosons at the tree level. However, even when the 
Higgs boson masses are set to be very large, the phase in the CKM matrix can be made 
finite and CP violating effects will not disappear unlike in other models of spontaneous CP 
violation where the CP violation disappear when Higgs boson masses become large. Another 
interesting feature of this model is that the FCNC Yukawa couplings are fixed in terms of 
the quark masses and CKM mixing angles, making phenomenological analysis much easier. 

We have studied some implications for meson-anti-meson mixing, including recent data 
on D — D mixing, and neutron electric dipole moment. We find that the neutral Higgs 
boson masses can be at the order of 100 GeV. The neutron EDM can be close to the present 
experimental upper bound. 
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